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Based on the t — t' — U — V model with proper chosen parameters for describing the cuprate 
superconductors, it is found that near the optimal doping at low temperature (T), only the pure 
d-wave superconductivity (dSC) prevails and the antiferromagnetic (AF) order is completely sup- 
pressed. At higher T, the AF order with stripe modulation and the accompanying charge order 
may emerge, and they could exist above the dSC transition temperature. We calculate the local 
differential tunnelling conductance (LDTC) from the local density of states (LDOS) and show that 
their energy variations are rather different from each other as T increases. Although the calculated 
modulation periodicity in the LDTC/LDOS and bias energy dependence of the Fourier amplitude of 
LDTC in the "pseudogap" region are in good agreement with the recent STM experiment [Vershinin 
et al., Science 303, f995 (2004)], we point out that some of the energy dependent features in the 
LDTC do not represent the intrinsic characteristics of the sample. 
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Recently, STM or the local differential tunnellingcon- 
ductance (LDTC) measurement by Vershinin et al. [jj] on 
slightly underdoped BSCCO indicated that the electronic 
states at low energies and at a temperature T higher than 
the superconducting transition temperature (T c ) in the 
pseudogap region exhibit an energy-independent spatial 
modulation which resulting to a checkerboard pattern 
with incommensurate periodicity 4.7a ±0.2 (a is the lat- 
tice constant). At very low temperature, however, no 
such pattern has been detected in agreement with 
previous measurements 0, 0. In addition, the Fourier 
amplitude of the LDTC at the modulation wave-vector 
which corresponds the periodicity 4.7a increases its mag- 
nitude first and then flattens out as the bias energy de- 
creases. How to understand these behaviors are out- 
standing questions which have not been addressed in the 
existing literatures. In the present paper we are trying 
to explain these issues by adopting the idea of the pos- 
sible coexistence of the d-wavc superconductivity (dSC) 
with antiferromagnetic (AF) order for cuprate supercon- 
ductors 0, 0, @j 1 an d to examine the formation of the 
AF order and the accompanying charge order at finite 
temperature. The phcnomcnological t — t' — U — V model 
will be applied to describe the cuprate superconductors. 
With proper chosen parameters, we show that at low 
temperature only dSC prevails in our system and the 
AF order is completely suppressed. The local density 
of states (LDOS) and LDTC images are featureless. At 
higher temperature, it is found that the AF order with 
stripe modulation, which is also referred to the spin den- 
sity wave (SDW), and the accompanying charge order 
or the charge density wave (CDW) may show up and 
they could even persist at temperatures above the BCS 
superconducting transition temperature T^ cs . In the 
presence of SDW, we show both of the LDOS/LDTC im- 
ages to have energy-independent stripe modulation with 



spacing of 5a/4a spreading over a 48a x 24a lattice. Ac- 
cording to the Fourier analysis of the LDOS images, an 
average periodicity 4.8a could be assigned for the stripe 
modulation. If both of the doubly degenerate x- and 
y- oriented stripes have the probability to appear in the 
time interval of the measurement or the proximity effect 
[U exists between neighboring domains with differently 
oriented stripe modulations, the combined LDOS images 
would have a checkerboard pattern of 4.8a x 4.8a struc- 
ture. All these features are consistent with the experi- 
ment of Vershinin et al. . In order to compare with the 
energy variation of the STM measurements at finite tem- 
perature (T), the LDTC is needed and it can be obtained 
from the LDOS by using the method of convolution, al- 
though both of the energy variations of LDTC and LDOS 
exhibit the "pseudogap" -like characteristics their be- 
haviors at higher temperature are quite different. While 
the size of the gap for quasiparticle excitations as a func- 
tion of T is measured by the separation between the co- 
herent peaks in the LDOS, the use of the LDTC or the 
experimental data to directly determined this quantity 
could be misleading. 

To model these observed phenomena, we employ an 
effective mean-field t — t' — U — V Hamiltonian by as- 
suming that the on-site repulsion U is responsible for the 
competing antifcrromagnctism and the nearest-neighbor 
attraction V causes the d-wave superconducting pairing 
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where is the hopping integral, /i is the chemical po- 
tential, and Ay = -|-( Ci T c ji — c U c jT) i s the spin-singlet 
d-wave bond order parameter. The Hamiltonian above 
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shall be diagonalizcd by using Bogoliubov-de Gennes' 
(BdG) equations, 
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where Hy^ = — fy + (U (n\ a ) — ju)<Jy. Here, we choose 
the nearest-neighbor hopping (ty) =4=1 and the next- 
nearest-neighbor hopping (fy) = t' = —0.25 to match 
the curvature of the Fermi surface for most cuprate su- 
perconductors . The exact diagonalization method to 
self-consistently solve BdG equations with the periodic 
boundary conditions is employed to get the N positive 
eigenvalues (E n ) with eigenvectors (u^, v^) and iV nega- 
tive eigenvalues (E n ) with eigenvectors (— The 
self-consistent conditions are 
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where u" = (—v^*,u^) and v™ = (ujj*,?;,™) are the row 
vectors, and f(E) = l/(e l3E + 1) is Fermi-Dirac dis- 
tribution function. Since the calculation is performed 
near the optimally doped regime, the filling factor, nf = 

Xd<r( c L CiCT ) /N x N y , is fixed at 0.85, i.e., the hole doping 
8 = 0.15. Each time when the on-site repulsion U or the 
temperature is varied, the chemical potential \i needs to 
be adjusted . 

In Fig. 1(a), we use U = 2.44 and V = 1.0 to approx- 
imately reproduce the phase diagram of Inaba et al. Q 
based on a mean field approach of the t — J model, there 
the contribution from the AF order was explicitly consid- 
ered. We shall apply this phase diagram to understand 
the experimental observations in Ref. 0] . It is important 
to notice that the curve (dashed line) separating the re- 
gion of the dSC and the coexistence region of dSC plus 
SDW in the phase diagram has a positive slop as the dop- 
ing ( 5) increases. The detailed structure of the SDW 
in the underdoped region is not the main focus of the 
present paper and will not be presented here. From the 
phase diagram, if <5 = 0.15 is chosen near the optimal dop- 
ing or in the slightly underdoped region, one can easily 
see that at low temperature T, our system is in the phase 
of dSC. When T increases, the SDW order emerges and 
coexists with the dSC. At very higher T, the dSC disap- 
pears and only SDW order survives. These temperature 
dependencies can be understood from Fig. 1(b) where 
the dSC order parameter at site i and the maximum value 
of the staggered magnetization Mi are plotted as func- 
tions of temperature T/t. From Fig. 1(b), it is straight- 
forward to show that only dSC exists and the AF order is 
completely suppressed at low T. As T > 0.06t, the stripe- 
modulated AF or the SDW order incipit, and both of 
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FIG. 1: (a) The phase diagram for U = 2.44. (c) The projec- 
tion of the staggered magnetization Mi at T = 0.84Tjv, where 
Mi = (— 1) (nif — nil). The size of the unit cell is N x x N y = 
48 x 24. (b) and (d) show the temperature dependence of the 
dSC (open circle) and the maximum value of Mi (solid circle) 
for U = 2.44 and U = 2.39, respectively. The value of the 
dSC order parameter Ai = ^(Ai+4 + Ai_£ — Ai+g — Ai_y) is 
measured in the unit of t. 



them could persist above the BCS transition temperature 
T C BCS = O.Ut. At T > Tf cs , the staggered magnetiza- 
tion decreases rapidly to zero at the Neel's temperature 
Tn = 0.155t. In Fig. 1(c) we show that the projection 
of the y-oriented stripe modulation in the AF order (or 
SDW) at T = 0.84T/V < T^ cs along x-axis, where the 
SDW has mixed periodicity 10a — 9a — 10a — 9a — 10a 
over a 48a x 24a lattice. The periodicity of the SDW 
seems not sensitive to T as long as the SDW order is in 
presence. In Fig. 1(d), we make similar plots as those 
in Fig. 1(b) but with a smaller U = 2.39. In this case 
the AF order is completely suppressed and our system 
is in the state of pure dSC at all temperatures. This 
implies that the phase diagram for U = 2.39 should be 
similar to that in Fig. 1(a) except the phase boundary 
(dashed line) is pushed toward the lower doped region. 
In both figures 1(b) and 1(d), the dSC order parameter 
as a function of T appears to have the BCS-likc behavior. 
Besides, the T ECS in the case for U = 2.39 [Fig. 1(d)] 
is slightly larger than the one for U = 2.44 [Fig. 1(b)]. 
This is because the appearance of SDW in Fig. 1(b) at 
higher T also suppresses the dSC. It needs to point out 
here that a phase diagram of the same model was previ- 
ously studied by Martin et al. using a very larger U 
on a much smaller lattice (17a x 10a). The curve sepa- 
rating the region of dSC and that of dSC plus SDW in 
their phase diagram has a negative slop. And that would 
yield the conclusion that if the system is in pure dSC 
state at low temperature, then it is always in dSC state 
at a temperature up to T < T^ cs and the SDW order 
never shows up. This is very different from the present 



3 



situation. In the following we show that the results in 
Figs. 1(b) and 1(c) can be applied to understand various 
features observed in the STM experiment 0. 

At finite temperature, what the STM measures is the 
LDTC which has the following definition 
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is the LDOS, and A is proportional to the square of the 
tunnelling matrix element. The summation in pi{E)t 
is averaged over M x x M y wavevectors in first Brillouin 
Zone. 

In Fig. 2, the normalized LDOS (a) and LDTC (b) 
at the sites with the maximum staggered magnetization 
(Mi) as functions of the bias energy E are presented from 
low to high temperatures. It is clear from Eq. (4) that 
the LDTC reduces to LDOS at low T. But at higher T, 
the behavior of the LDTC is rather different from that of 
the LDOS. For example, the characteristic temperature 
rpbuik ^ w ] :iere coherent peaks of the superconductivity 
begin to flatten out, in the LDTC is considerably lower 
than that in the LDOS. One usually define the gap of 
the quasiparticles as the separation between the coher- 
ent peaks or the width of the dip or depression in the 
LDOS. While the gap in the LDOS [Fig. 2(a)] appears 
roughly to be a constant from low T up to T/v, the "gap" 
in the LDTC [Fig. 2(b)] increases progressively from low 
T to high T. The increment in the magnitude of "gap" 
as T raises has been indeed observed by experimental 
measurements 0, E3> E3 • It should be emphasized here 
that the "gap" in the LDTC or observed directly from 
the STM experiments is the result of the convolution in 
Eq. (4) and may not be the real gap of the system. In the 
temperature range of T^ ulk < T < Tn as shown in Fig. 
2, the system appears to be in the "pseudogap" region 
according to the characterization from STM experiments 
0, Q. Even though the effect due to the phase fluc- 
tuations on the dSC order parameter 0, 0] has not 
been taken into account. In this region, the Fermi sur- 
face in our theory is everywhere gaped while the ARPES 
experiment [l^ indicates that the gaps occur only near 
(±7r,0) and (0, ±7r). This difficulty so far has not been 
understood. 

It is well established that associated with the SDW 
shown in Fig. 1 (c), there are stripe- modulations in the 
dSC order and the charge density. The stripe modulated 
charge order can also be refereed to the charge density 
wave (CDW). In the top two graphs of the left panel in 




FIG. 2: Temperature dependence of the normalized 
LDOS (pi(E) T /pi(E) TN ) and the normalized LDTC 
{Gi{E)t/Gi(E) Tn ), here T N = 0.155*. The LDOS and 
LDTC at Tat as a function of E is represented by the dashed 
line. This representative set of spectra was shifted vertically 
for clarity. The wavevectors in first Brillouin Zone are 
M x x My = 24 x 12. 



Fig. 3, the projections of the y-oriented stripe modula- 
tions in the dSC order and the charge density are plotted 
along the z-axis at T = 0.84Tjv(< T® cs ). Here the spa- 
tial distributions of the dSC and the charge order are only 
slightly modulated by the stripe structure, and the dSC 
and the SDW orders coexist in real space. This feature is 
very different from the case in Ref. 13, where the rfSC or- 
der is practically suppressed to zero in the spatial regions 
where the AF order is in presence. The stripe modula- 
tion also appears in the LDTC (or LDOS) as shown in 
the bottom graph of the left panel. The stripe modu- 
lations displaced in the left panel have the same mixed 
periodicity 5a — 5a — 4a — 5a — 5a over a 24a lattice 
along a;-axis. This indicates that the system is trying to 
establish a periodicity incommensurate with the under- 
line lattice, but fails to do so because the calculation is 
performed in a discrete and finite lattice, not in a con- 
tinuum. It is also straightforward to show that the same 
mixed periodicity still remains in the CDW and LDTC 
(or LDOS) even T is in the region of Tf cs <T <T N , 
where the superconducting order parameter vanishes and 
the SDW order is still in presence. As it will be shown 
below that an average periodicity 4.8a for the stripe mod- 
ulation could be assigned in the present case. The x- and 
y- oriented stripe modulations are degenerate in energy, 
and it is thus possible for the both x- and y- oriented 
stripe-modulations to show up either in the time duration 
when the experiment is performed or due to the proxim- 
ity effect between neighboring domains with stripes of 
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FIG. 3: Left panel: From top to bottom, the projection of 
the dSC order, charge density, and LDTC at T = 0.84Tjv. 
Right panel: The energy evolution of the normalized Fourier 
amplitude of the LDOS (a) and LDTC (b) images for T = 
0.94Tjv (solid circle) and T = 0.86Tjv (open circle) at q x = 
0.417(7r/o). The solid curve represents the normalized LDOS 
(a) and LDTC (b) as a function of E at T = 0.84T N . The 
wavevectors which have been used to calculate the LDOS and 
LDTC images are M x x M y = 6 x 6 in first Brillouin Zone. 



different orientations. As a result, checkerboard pattern 
could be observed. 

Furthermore, in order to get a detailed understanding 
of the periodicity of the modulation in the LDOS (or 
LDTC) image, we perform the Fourier transform of the 
LDOS (or LDTC) image in the temperature range where 
the SDW order is in existence. 
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^exp(zq-ri) ■ pi(E) T . (6) 



When the above quantity is plotted against q x , a sharp 
peak occurs around q x = 0A17(n/a) = 27r/4.8a, indicat- 
ing an average "periodicity" 4.8a in real space. In Fig. 
3(a) and 3(b), the solid curves at the bottom represent, 
respectively, the normalized LDOS and the normalized 
LDTC as functions of the bias energy at T = 0.84Tjv. 
The curves at the top made of solid/open circles show 
the bias energy dependencies of the Fourier amplitudes 
of the LDOS and LDTC images at q x = 27r/4.8a for 
two different temperatures normalized by their values 
at E = 0.6i. The Fourier amplitude of the LDOS first 
reaches a peak and then drops somewhat rapidly as the 
bias decreases to zero. Near the zero bias, the result 
at T = 0.8AT N {< T C BCS ) dips more than the result at 
T = 0.94T W (> T C BCS ) in the pseudogap region. On the 
other hand, the Fourier amplitude of the LDTC as func- 
tions of the bias energy at these two temperatures differ 
very little, and they drops only slightly after reaching a 
broader maximum from higher bias. This feature in the 



LDTC seems to be in better agreement with the STM 
measurements Q as compared with that of the LDOS. 

In conclusion, the temperature dependencies of the 
LDTC and LDOS in a cuprate superconductor with the 
competing AF order have been investigated in the present 
paper. According to our calculations based on the phase 
diagram in Fig. 1(a) near the optimal doping, there is 
no signature of the charge order in the LDOS and the 
LDTC at low T. When the SDW order appears, both 
of the LDOS and the LDTC exhibit the same CDW like 
modulation with an average periodicity 4.8a at T below 
and above T BCS . We also calculate the Fourier ampli- 
tude of LDOS and LDTC at q = (27r/4.8a,0) and their 
bias energy dependencies in the "pseudogap" region as 
shown in Figs. 3(a) and 3(b), respectively. All the fea- 
tures in the LDTC are in good agreement with those ob- 
served by the STM experiments [l| . Finally we point out 
that the " gap" of the quasiparticles at higher T obtained 
directly from the LDTC or the STM experiments does 
not correspond to the real gap of the system in cuprate 
superconductors. 
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